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ABSTRACT 

This paper uses a single unifying principle (based upon the nonlinear “momentum-flux” 
effects produced when different components of a motion transport different components of 
its momentum) to give a broad scientific background to several aspects of the interaction 
between airflows and atmospheric sound. First, it treats the generation of sound by airflows 
of many different types. These include, for example, jet-like flows involving converted turbu- 
lent motions with the resulting aeroacoustic radiation sensitively dependent on the Mach 
number of convection - and they include, as an extreme case, the supersonic “boom” (shock 
waves generated by a supersonically converted flow pattern). Next, the paper analyses sound 
propagation through nonuniformly moving airflows, and quantifies the exchange of energy 
between flow and sound; while, finally, it turns to problems of how sound waves “on their 
own” may generate the airflows known as acoustic streaming. 


'This research was supported by the National Aeronautics and Space Administration under NASA (ton- 
tract No. NAS1-18605 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23681-0001. Research 
was also supported by the Leverhulme Trust. 
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1. Broad Overview 


This general introductory paper is devoted to Interactions of Sound with Air, including 
transmission through the atmosphere and 

both generation of sound by 
and propagation of sound in 

(e.g,, manmade flows - around aircraft or air machinery - or natural winds) as affected by 
the airs boundaries and atmospheric composition; with (conversely) generation of airflows 
by sound (acoustic streaming). 

From linear acoustics I utilize the properties of the wave equation, including 

( 1 ) the short-wavelength ray-acoustics approximation (Lighthill, 1978b, hereinafter de- 
noted WF, pp. 67-) and 

(ii) the theory of multipole sources (WF, pp. 31-) - with the long- wavelength “compact- 
source approximation (source region of size f with u>(/c small, where o> = radian 
frequency, radiates like a concentrated source); 

while from nonlir-ar acoustics I use (WF, pp. 150 ) the physics of waveform shearing and 
shock formation. 

Techniques special to aeroacoustics and atmospheric sound are centered on the momen- 
tum equation for air. Its difference from a wave-equation approximation include 

A. Linear effects, of gravity acting on air stratified as meteorologists observe; effects which 
allow independent propagation (WF, pp. 292-) of “internal” gravity waves and of 
sound, except at wavelengths of many kilometers when the atmosphere becomes a 
waveguide (WF, pp. 425-) for global propagation of interactive acoustic-gravity waves; 

and (still more importantly) include 

B. Nonlinear effects, of the momentum flux pu iU] ; i.e. the flux - rate of transport across 
unit area - of any pu, momentum component by any u } velocity component. This 
term, neglected in linear acoustics, acts like a stress (i.e. force per unit area - since 
rate of change of momentum is force). In particular, 

(i) an airflow’s momentum flux pu,u ; generates sound like a distribution of (time- 
varying) imposed stresses; thus not only do forces between the airflow and its 
boundary radiate sound as distributed dipoles, (-> F ~ ~ + F) but also such 
stresses (acting on fluid elements with equal and opposite dipole-like forces) ra- 
diate (WF, pp. 63-) as distributed quadrupoles; (-♦«— ~ — (- + _) (Lighthill, 
1952 and Lighthill, 1962) 
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(ii) the mean momentum flux < piiiUj > in any sound waves propagating through 
a sheared flow (with shear dVi/dxj) is a stress on that flow (Lighthill, 1972 and 
WF, pp. 329-), and the consequent energy exchange (from sound to flow when 
positive, vice versa when negative) is 

< puiUj > dVi/dxj] ( 1 ) 

(iii) even without any pre-existing flow, energy-flux attenuation in a sound wave al- 
lows streaming to be generated (WF, pp. 337-) by unbalanced stresses due to 
a corresponding attenuation in acoustic momentum flux - essentially, then, as 
acoustic energy flux is dissipated into heat, any associated acoustic momentum 
flux is transformed into a mean motion (Lighthill, 1978a and WF). 

And another (less crucial) momentum-equation/wave-equation difference is 

C. Nonlinear deviation of pressure excess p - p 0 from a constant multiple, c 2 0 {p - p 0 ), of 

density excess. 

(i) For sound generation by airflows, this adds an isentropic term to the quadrupole 
strength per unit volume (WF, pp. 60-) 

Tij = P U > U J + [(? - Po) - clip - Po)] Sij, (2) 

the last term being considered important mainly for flows at above-ambient tem- 
peratures (Lighthill, 1952 and 

Lighthill, 1963); 

(ii) for propagation of sound with energy density E through flows of air with adiabatic 
index 7, the mean deviation is about (7 - 1 )E, and the total radiation stress 
(Bietherton and Garrett, 1968) 

< puiUj > +^(7 - \)ESij (3) 

adds an isotropic pressure excess to the mean momentum flux (although the 
energy exchange (1) is unchanged in typical cases with dl’i/dx, essentially zero). 

(And we note that the very special case of sound waves interacting on themselves [in other 
words, nonlinear acoustics] may be interpreted (WF, p. 148) as a combined operation of the 
“self-convection effect” B and the (smaller) “sound-speed deviation” C.) 
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when flow 
instabilities 
lead to 


regular flow 
oscillations; 
extremely irregular 
turbulence. 


2. Compact Source Regions 

2.1. Sound generation by low-Mach-number airflows 

The main nondimensional parameters governing airflows of characteristic speed U and 
length-scale L are Mach number M = U/c (where, in aeroacoustics, c is taken as the sound 
speed in the atmosphere into which sound radiates) and Reynolds number R = UL/u where 
kinematic viscosity. Low-Mach-number airflows are compact sources of sound, with 
frequencies 

narrow-banded 
at moderate R 
broad- banded 
at high R 

Since, in either case, a typical frequency u> scales as U/L (Strouhal scaling), the compactness 
condition uL/c small is satisfied if M = U/c is small (Lighthill, 1962). 

A solid body which, because of flow instability, is subjected to a fluctuating aerodynamic 
force F scaling as pU 2 L 2 (at frequencies scaling as U/L), radiates as an acoustic dipole of 
strength F, with mean radiated power < F 2 > /127rp.: 3 . 

This acoustic power scales as pU 6 L 2 /c 3 (a sixth-power dependence on flow speed). There- 
fore acoustic efficiency, defined as the ratio of acoustic power to a rate of delivery (scaling 
as pU 3 L 2 ) of energy to the flow, scales as (U/c) 3 = M 3 . 

(Exceptions to compactness include bodies of high aspect-ratio; thus, a long wire in 
a wind [where the scale L determining frequency is its diameter] radiates as a lengthwise 
distribution of dipoles.) 

Away from any solid hody a compact flow (oscillating or turbulent, with frequencies 
scaling as U/L) leads to quadrupole radiation (see B(i) above) with total quadrupole strength 
scaling as pU 2 L 3 . Acoustic power then scales as pU s L 2 /c 5 : an eighth-power dependence 
(Lighthill, 19o2 and Lighthill, 1962) on flow speed. In this case acoustic efficiency (see 
above) scales as (U/c) 5 — A/ 5 . 

Such quadrupole radiation, though often important, may become negligible near a solid 
body when dipole radiation due to fluctuating body force (with its sixth-power dependence) 
is also present (Curie, 1955 and Ffowcs Williams and Hawkings, 1969). 

Near not necessarily compact bodies a more refined calculation - using Green’s functions 
not for free space but for internally bounded space leads in general to the same conclu- 
sion: that quadrupole radiation with its eighth-power dependence is negligible alongside the 
sixth-power dependence of dipole radiation due to fluctuating body forces; but important 
exceptions to this rule include sharp-edged bodies, where features of the relevant Green’s 
function imply a fifth-power dependence on flow speed of acoustic radiation from turbulence 
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(Ffowos Williams and Hall, 1970; Crighton and Leppington, 1970; Crighton and Leppington, 
1971; and Crighton, 1981). 

2.2. Sound generation by turbulence at not so low Mach number 

The chaotic character of turbulent flow fields implies that velocity fluctuations at points 
P and Q , although they are well correlated when P and Q are very close, become almost 
uncorrelated when P and Q are not close to one another. 

Reminder: statisticians define correlation coefficient C for the velocities Up and Uq as 
C =< vp • vq > / < Vp >’/ 2 < \q in terms of the deviations, \p = Up— < Up > 
and Xq = Uq— < Uq >, from their means. When two uncorrelated quantities are combined, 
their mean square deviations are added up: 

< Vp + Vq > 2 = < Vp > +2 < Vp • Vq > + < Vq > 

= < Vp > + < Vq > if C = 0. 

Theories of turbulence define a correlation length (, with 

. | well correlated 1 . | 1 1 

and u « \ uncorrelated / “ ciose lo { 0 / " ~ 

when PQ is substantially 

Roughly speaking, different regions of size ( (“eddies”) generate sound independently, and 
the mean square radiated noise is the sum of the mean square outputs from all the regions 
(Lighthill, 1954). Typical frequencies in the turbulence are of order = v/( , where v 
is a typical root mean square velocity deviation < v 2 > *' 2 , so that for each region the 
compactness condition wf/c small is satisfied if v/c is small. Compactness, then, requires 
only that a r.in.s. velocity deviation v (rather than a characteristic mean velocity U) be 
small compared with c - which is less of a restriction on M — V/c and can be satisfied at 
“not so low” Mach number (Lighthill, 1963). 

3. Doppler Effect 

How is the radiation from such “eddies” modified by the fact that they are being converted 
at “not so low” Mach number? The expression Doppler effect, covering all aspects of how 
the movement of sources of sound alters their radiation patterns, comprises (i) frequency 
changes (WF), (ii) volume changes (Lighthill, 1952 and Lighthill, 1962), (iii) compactness 
changes (Lighthill, 1963 and Ffowcs Williams, 1963). 
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3.1. Frequency changes 


When a source of sound at frequency u> approaches an observer at velocity w , then in a 
single period 1 = 2it/u sound emitted at the beginning travels a distance cT while at the 
end of the period sound is being emitted from a source that is closer by a distance wT. The 
wavelength A (distance between crests) is reduced to 


A = cT — xvT = ‘27r(c — w)/u> 


(4) 


(Figure 1) and the frequency heard by the observer (2 tt divided by the time A/c between 
arrival of crests) is increased to the Doppler-shifted value (WF) 


U! r = 


LO 


- («V0 


the relative frequency 


(5) 


that results from relative motion between source and observer. For an observer located on a 
line making an angle 6 with a source’s direction of motion at speed V, the source’s velocity of 
approach towards the observer is w = V cos 6 (Figure 1 ) and the relative frequency becomes 


U 


LJr — 


{ augmented 
diminished 


when 6 is an 


1 ~ ( V/c ) cos 6 

Such Doppler shifts in frequency are familiar everyday experiences. 


{ acute ) 
obtuse J 


angle. 


( 6 ) 


3.2. Volume changes 

When an observer is approached at velocity w by a source whose dimension (in the 
direction of the observer) is f, sounds arriving simultaneously (Figure 2) from the source’s 

{ near } sides have been emitted ( YakT } b y a t (say). 

In the time t for sound from the far side to reach the observer, after travelling a distance 
c/, the relative distance of the near side in the direction of the observer was increased from 

t to f + ujr before it emitted sound which then travelled a distance c(t - r). Both sounds 
arrive simultaneously if 

ft = e u;t + c(«: - r), giving r = 

c — w 

( ^ 

and ( + wt - — = (jjJu. 

1 - (w/c) 

The source’s effective volume during emission is inereased, then, by the Doppler factor u> r /u> 
(since dimension in the direction of the observer is so increased whilst other dimensions are 
unaltered: Lighthill, 19. r )2 and Lighthill, 1962). 


If turbulent “eddies” are effectively being convected, relative to the air into which they 
are radiating, at velocity V, then Equation (6) gives, for radiation at angle 0 , the Doppler 
factor u> r /u> which modifies both the frequencies at which they radiate and the effective 
volume occupied by a radiating eddy. 

But Equation (2) specifies the quadrupole strength T t; per unit volume for such an eddy. 
Without convection the pattern of acoustic intensity around a compact eddy of volume P 
and quadrupole strength PT t] would be 

< (? 3 f i} XiXjr~ 2 )' 2 > /167tV/)oc 5 ; (8) 

and, since different eddies of volume P radiate independently, we ran simply add up mean 
squares in the corresponding expressions for their far-held intensities. This gives 

P < {fijxmr- 2 ) 2 > /167r 2 rV 5 (9) 

as the intensity pattern radiated by unit volume of turbulence. The Doppler effect modifies 
this, when the compactness condition is satisfied, by five factors u> r /u; (one for the change 
in source volume P and four for the frequency change as it affects the mean square of a 
multiple of the second time-derivative of T l} ) and this intensity modification by a factor 

[1 - (Vyc)cos 0] -5 (10) 

brings about an important preference for forward emission (Ffowcs Williams, 1963). 

3.3. Compactness changes 

As ( V/c ) increases, however, the Doppler effect tends to degrade the compactness of 
aeroacoustic sources in relation to forward emission. Not only does u>(/c increase in propor- 
tion to Mach number, but an even greater value is taken by u> r (/c, the ratio which must be 
small if convected sources are to be compact. A restriction on the extent (10) of intensity 
enhancement for forward emission as V/c increases is placed by these tendencies (Lighthill, 
1963; Ffowcs Williams, 1963; and Dowling et al., 1978). 

They can develop, indeed, to a point where the compact-source approximation may 
appropriately be replaced by its opposite extreme: the ray-acoustics approximation (WF). 
Thus, for supersonic source convection (V/c > 1), the relative frequency (6) becomes infinite 
in 

the Macli direction 6 = cos -1 (c/l/), (11) 

and radiation from the source proceeds (Figure 3) along rays emitted at this angle (Ffowcs 
Williams and Maidanik, 1965). 
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Explanatory note: the source’s velocity of approach w towards an observer positioned 
at an angle (11) to its direction of motion is the sound speed c; thus, not only is the generated 
wavelength (4) reduced indefinitely (the ray-acoustics limit) but, essentially, different parts 
of a signal are observed simultaneously: the condition of stationary phase satisfied on rays 
(VVF). Sounds emitted (WF, p. 196) by a source approaching an observer at a speed w 
exceeding c are heard by him in reverse order (“pap pep pip pop pup” becomes “pup pop 
pip pep pap”!) but when w = c all the sound (vowels and consonants!) are heard together 
as one single “boom.” 

Further note: the influences placing a limit on the signal propagated along rays may 
include the duration 8 of well-correlated emission from turbulent “eddies”; and, also, may 
include nonlinear effects (see §4.2. Supersonic booms). 


3.4. Uniformly val ; d Doppler-effect approximations 


Just as a correlation length ( for turbulence was specified in §2.2, so a correlation duration 
8 can be characterized by the requirement that moving eddies have { We ^ corr ^ a ^ 1 


J M { uncorrelated J 

velocities at times differing by substantially j < ^ |. Combined use of correlation length ( 

and duration 8 affords an approximation to the radiation pattern from convected “eddies” 
that has some value at all Mach numbers, spanning the areas of applicability of the compact- 
source and ray-acoustics approximations. 

Figure 4 uses space-time diagrams where the space-coordinate (abscissa) is distance in the 
direction of the observer. Diagram (a) for unconverted “eddies” approximates the region of 
good correlation as an ellipse with axes ( (in the space direction) and 8 (in the time direction). 
Diagram (b) shows such a region for convected “eddies” whose velocity of approach towards 
the observer is w\ thus, it is Diagram (a) sheared by distance w per unit time. 

Signals from far points F and near points N , in either case, reach the observer simulta- 
neously - as do signals from other points on the line FN - if this line slopes by distance c 
(the sound speed) per unit time. 

Compact-source case (i) with w/c small: the space component of FN in Diagram (b) is 
f[l — (u;/c)] _l , just as in Equation (7) for normal Doppler effect (neglecting finite 8). 

Ray-acoustics case (ii) with w/c = 1: the space component of FN is c8. 

Intermediate case (iii) with w/c “moderately” < 1: the space component of FN is ( 
multiplied by an enhancement factor 


[(l-te/c^-Mf/c*) 2 ] 


,AVl~>/2 


which represents the effective augmentation of saurce volume due to convection (Ffowcs 


0 . 
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Williams, 1963). 

This enhancement factor ( 12) is applied not only to the volume term f 3 in the quadrupole 
field (9) but also twice to each of the pair of twice-differentiated terms inside the mean square; 
essentially, because time-differentiations in quadrupole fields arise (WF) from differences 
in the time of emission by different parts of the quadrupole source region (and the time 
component of FN in Diagram (b) is simply the space component divided by c). As before, 
then, five separate factors ( 12) enhance the intensity field; and, with w replaced by V cos 0, 
expression (10) for the overall intensity modification factor is replaced by 

{[1 - (V/c)cos$} 2 + (e/c6) 2 }~ &/2 . (13) 

This modification factor (13) affords us an improved description of the influence of Doppler 
effect not only on the preference for forward emission but also on the overall acoustic power 
output from convected turbulence (Light-hill, 1963 and Ffowcs Williams, 1963). For example, 
Diagram (c) gives (plain line) a log-log plot of the average (spherical mean) of (13) as a 
function of V/c on the reasonable assumption that i — 0.61/ S. As V/c increases this average 
modification factor rises a little at first, but falls drastically like 5(1'/ c) -5 for V/c significantly 
> 1 . 

Now low- Mach- number turbulence away from solid boundaries (§2.1) should radiate 
sound with an acoustic efficiency scaling as (U/c.) 4 5 where U is a characteristic velocity in the 
flow. With V taken as that characteristic velocity (although in a jet a typical velocity V of 
eddy convection would be between 0.5 and 0.6 times the jet exit speed), the modification of 
(say) an acoustic efficiency of lO'^V'/c) 5 for low Mach number by the average modification 
factor would cause acoustic efficiency to follow the broken-line curve in Diagram (c), tending 
asymptotically to a constant value, 0.005, (aeroacoustic saturation) at high Mach number. 
Such a tendency is often observed for sound radiation from “properly expanded” supersonic 
jets (see below). 

4. Introduction to Aircraft Noise 
4.1. Aero-engine and airframe no ise 

How are aeroacoustic principles applied to practical problems - such as those of studying 
aircraft noise with a view to its reduction (Orighton, 1975; Goldstein, 1976; and Goldstein, 
1984)? 

In any analysis of the generation of sound by airflows, we may need first of all to ask 
whether the geometry of the problem has features that tend to promote resonance. For 
example, a long wire in a wind (§2.1) generates most sound when vortex-shedding frequencies 
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are fairly close to 
and so can “lock on” to 


the wire’s lowest natural frequency of vibration; giving good 


correlation of sideforces, and so also of dipole strengths, all along the wire. 

Again, a jet emerging from a thin slit may interact with a downstream edge (parallel 
to the slit) in a resonant way (Curie, 1953 and Powell, 1953b); with very small directional 
disturbances at the jet orifice being amplified by flow instability as they move downstream 
to the edge, where they produce angle-of-attack variations. Dipole fields associated with the 
resulting sideforces can at particular frequencies reneA' the directional disturbances at the 
orifice with the right phase to produce a resonant oscillation. Some musical wind instruments 
utilize such jet-edge resonances, reinforct 1 by coincidence with standing-wave resonances in 
an adjacent pipe. 

But in the absence of such resonances (leading to enhanced acoustic generation at fairly 
well defined frequencies) airflows tend to generate acoustic “noise” whose reaction on the 
flow instability phenomena themselves is negligible. 

Resonances analogous to the above which need to be avoided in aircraft design include, 
for example, 

(a) panel flutter, generated at a characteristic frequency as an unstable vibration of a 
structural panel in the presence of an adjacent airflow (Dowell, 1975); 

(b) screeching of supersonic jets from nozzles which, instead of being “properly expanded” 
so that an essentially parallel jet emerges, produce a jet in an initially non-parallel form 
followed by shock waves in the well known recurrent “diamond” shock-cell pattern; the 
first of these, replacing the edge in the above description, can through a similar feedback 
of disturbances to the jet orifice generate a powerful resonant oscillation (Powell, 1953a; 
Powell, 1953b; and Howe and Ffowcs Williams, 1978). 

And undesirable resonances may also be associated with aeroengine combustion processes 
(Candel and Poinsot, 1988). But we turn now to the aircraft noise of a broad-banded 
character that remains even when resonances have been avoided. 

Then aero-engine jet noise proper ( Lighthi II , 1963) (thet is, the part unrelated to any 
interaction of jet turbulence with solid boundaries) tends to follow a broad trend similar to 
that in Figure 4; where, however, because the eddy convection velocity V is between 0.5 
and 0.6 times the jet exit speed U, the acoustic efficiency makes a transition between a 
value of around 1 0~ 4 A/ 5 in order-of-magnitude terms for subsonic values of M = U/c and 
an asymptotically constant value of 10~ 2 or a little less for M exceeding about 2. 

The above tendency for M < 1 implies that noise emission from jet engines may be 
greatly diminished if a given engine power can be achieved with a substantially lower jet 
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exit speed, requiring of course a correspondingly larger jet diameter, L. Furthermore, with 

acoustic power output scaling as V (§2- ' ) lhr " sl “ ^ ™"“ 0B ” 

given thrust can be greatly reduced if V can be decreased and L increased by compara , e 

factors. . . . 

Trends (along these lines) in aero-engine design towards large turbofan engines with 

higher and higher bypass ratios, generating very wide jets at relatively modest mean Mach 

numbers, have massively contributed to iet noise suppression (whilst also winning advantages 

of reduced fuel consumption). On the other hand, such successes in suppressing jet noise 

proper (originally, the main component of noise from jet aircraft) led to needs for a Cheated 

focusing of attention upon parallel reductions of other aircraft-noise sources (Criglitoii, 19,2): 

(a) those associated with the interaction of jet turbulence with solid boundaries where 
sharp-edged boundaries (§2.1) pose a particular threat, 

(b) fan noise emerging from the front of the engine and turbine noise from the rear: 

(c) airframe noise including acoustic radiation from bonndary-laver turbulence and from 
interaction of that turbulence with aerodynamic surfaces for control purposes or lift 

enhancement. 

q 0 .>-,c key areas of modern research on aero engine and airframe noise are: 
for jet noise, techniques for relating acoustic output to vorticity distributions (Powell 
1964 and Mohring, 1978), and to any coherent structures (Kilmer, 1964 and Flours W illiams 
and Kempton, 1978). in jet turbulence; and for taking into account (cf § ,.4) propagation 

through the sheared flow in a wide jet (Phillips, 1960 and Maui, 1976): 

for noise from fans and propellers, mathematically sophisticated ways of reliable esti- 
mating the extent of cancellation of dipole radiation from different parts of a rotating blade 
system (alongside a good independent estimate of quadrupole radiation: Parry and ngliton. 

1991); , . , 

for airframe noise, a recognition (Powell, 1960 and (VigliUm, 1984) that massive <a„« - 

latinos act to minimize noise radiation from boundary-layer tiirbuh n<< on a flat sml 

uniform compliance - and, therefore, that avoidance of sharp no, .uniformities airframe 

skin compliance may promote noise reduction. 

4.2. Su personic booms 

1„ addition to aeroengine anil airframe noise, any aircraft flymg at a siipersonii spied 
v emits a concentrated •hoo.n’-like noise along rays (Figure ») in the Marl, direction ( ). 

I sketch the theory of supersonic booms with the atmosphere approximated as isotherma 
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{so that the undisturbed sound speed *akes a constant value c even though the undisturbed 
density p varies with altitude): a case permitting quite a simple extension of the nonlinear 
analysis of waveform shearing and shock formation ( WF and Whitham, 1956). Then the rays 
continue as straight lines at the Mach angle for reasons summarized in the explanatory note 
below expression (11). (Actually, the slight refraction of rays by temperature stratification 

m the atlT1 °sphere, when taken into account in a generalized version of the theory, produces 
only somewhat minor modifications of the results.) 

As such straight rays stretch out from a straight flight path along cones with semi-angle 
(11), any narrow tube of rays has its cross-sectional area A increasing in proportion to 
distance r along the tube (WF and Lighthill, 1956). On linear theory (WF), acoustic energy 
flux u 2 pcA is propagated unchanged along such a ray tube (so that u(pr) 1 / 2 is unchanged) 
where u is air velocity along it. On nonlinear theory, u(pr) 1/2 is propagated unchanged but 
at a signal speed altered to 

„ , 1+1 

+ ~J~ u (14) 

by self-convection and excess-wavespeed effects. 

This property can be described (WF) by an equation 


c 


7 + 1 u \ d d 
2 c*7 dt + dr 


u{pr) x/2 = 0 


(15) 


(WF, pp. 187-) where the quantity in braces is the altered value of the reciprocal of the 
signal speed (14). Now a simple transformation of variables 


x \ — r ct, t] — (pr) x/2 dr, u, = - -{pr) x/2 (16) 

is found to convert Equation (15) into the familiar form 

du\ dv.i 

aTV + u, a77 = 0 < l7 > 

which describes the waveform shearing at a uniform rate that is associated with shock for* 
mation and propagation in nonlinear plane-wave acoustics (Figure 5) 

From amongst this equation’s physically relevant solutions - namely, those with area- 
conserving discontinuities (representing shocks) - the famous /V- wave solution is the one 
produced by an initial signal (such as an aircraft’s passage through the air) that is first 
compressive and then expansive. The rules (WF) governing A -wave solutions of Equation 
(17) are that the discontinuity Au, at each shock fal's off like t^ ri while the space (change 
Ar, in x,) between shocks increases like t\ r \ These rules for the transformed variables (16) 
have the following consequences for the true physical variables: at a large distance r from 


the flight path the velocity change A u at each shock and the time interval At between the 
two shocks vary as 


Au 


(pi') j ( pr ) 1/,2 <fr and At ss J (pv) 


1/2 


(18) 


On horizontal rays (at the level where the aircraft is flying), p is independent of r and the 
Equations (18) take the greatly simplified form 


Au r 3 ^ 4 and At ~ r 1 ^ 4 


(19) 


appropriate to conical N - waves in a homogeneous atmosphere (which supersonically con- 
vected eddies may also generate). Actually, the rules (19) apply also to the propagation 
of cylindrical blast waves generated by an exploding wire; since, here also, ray tube areas 
increase in proportion to r. 

On downward pointing rays in an isothermal atmosphere p increases exponentially in 
such a way that the time interval At between shocks approaches the constant value obtained 
in (18) by making the integral’s upper limit infinite (WF and Lighthill, 1956). On the 
other hand the shock strength (proportional to the velocity change Au) includes the factor 
(pr) -1 / 2 where the large increase in p fiom the flight path to the ground (as well as in 
r) enormously attenuates the supersonic boom. Below Concorde cruising at Mach 2, for 
example, an observer on the ground hears two clear shocks with an interval of around 0.5s 
between them, and yet with strengths A p/p only about 0.001. 


5. Propagation of Sound through Steady Mean Flows 
5.1. Adaptations of ray acoustics 

Useful information on sound propagation (including sound of aeroacoustic origin) through 
steady mean flows (Blokhintsev, 1956 and Lighthill, 1972) can be obtained by adaptations 
of the ray-acoustics approximation. I sketch these here before, first, applying them (in §5.3 
below) to propagation through sheared stratified winds and, secondly, giving indications 
of how effects of such parallel mean flows are modified at wavelengths too large for the 
applicability of ray acoustics. 

Sound propagation through a steady airflow represents an autonomous mechanical sys- 
tem: one governed by laws that do not change with time. Then small disturbances can be 
Fourier-analyzed in the knowledge that propagation of signals with different frequencies ui 
must proceed without exchange of energy between them. 
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General Theory of Such Systems (WF, pp. 317-) 

Such disturbances of frequency u> involve pressure changes in the form Pcosa where P 
varies with position and the phase a is a function of position and time satisfying 

^ = w : the frequency; and - ^ = k, : the wavenumber, (20) 

a vector with its direction normal to crests and its magnitude 2ir divided by a local wave- 
length. 

In ray theory for any wave system (WF and Lighthill, 1972), we assume that the wave- 
length is small enough (compared with distances over which the medium - and its motion, 
if any - change significantly) for a well-defined relationship 

u> = n(k„Xi) (21) 

to link frequency with wavenumber at each position. Equations (20) and (21) require that 

_dk 1 _ d 2 a _ du _ dQ / d 2 o \ dQ_ _ dQ dk, dQ 

dt dxjdt dx, ~ dk, \ dxidx") + dx, ~ dk,~d7, + dx,' ^ 22) 

yielding the basic law (in Hamiltonian form) for any wave system: 

.. , . dx , dQ dk, dQ 

on rays satisfying — - = wavenumbers vary as — - = ; (23) 

dt ok x dt dxj 

equations easy to solve numerically for given initial position and wavenumber. However, the 
variations (23) of wavenumber ( “refraction” ) produce no change of frequency along rays: 


du _ dQ dk x dQ dx{ dQ / dQ 
~ “ ~~ + “ “ = dk, \ Ihi 


dt dk, dt dx, dt dk 


8Q dQ 
^ Oxi dk, 


so that rays are paths of propagation of the excess energy, at each frequency, associated with 
the waves’ presence. 

For sound waves we write k as the magnitude of t he wavenumber vector, expecting that 
at any point the value of the relative frequency in a frame of reference moving at the local 
steady flow velocity uj, will be cjk (the local sound speed times k)\ this implies (WF and 
Lighthill, 1972) that 

da da 

Ur ~ ~dt + Uj, d7, = u ~ U/,k " giving u==UJr + n !' k ' - c l k + «/ k > (25) 

as the acoustic form of the relationship (21). 

[Note: this rule (25) for relative frequency agrees with the Doppler rule (6), since the 
velocity of a source of frequency u> relative to stationary fluid into which it radiates is minus 
the velocity of the fluid relative to a frame in which the acoustic frequency is u;.] 
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Use of this form (25) of the relationship (21) in the basic law (23) tells us that 
dkj _ dcj dufi dx i- 

dt - - *• ^7 on ^ with sr = c Jj+ «/« (26) 

»l.ere the las. terms in these equations represent adaptations of ray acoustics associated 

" the mean f ’ ow - F °' the velocity of propagation along rays is the vector sum 

Of the mean flow velocity Ujl with a wave velocity of magnitude c, and direction normal to 
crests. 


5 Lner gy exchange between sound waves and mean flow 

The excess energy (say, E per unit volume) associated with the presence of sound waves 
is Propagated along such rays; in particular, if attenuation of sound energy is negligible, then 

flux of excess energy along a ray tube = constant. (27) 

Note; this excess energy density E is by no means identical with the sound waves’ energy 
density 6J 

E ’ =< 2^“"”" > + < iWP, >= c U'i' < !>]> (28) 

(Where the subscript s identifies changes due to the sound waves and the equality of the 
kmet.c and potential energies makes E, simply twice the latter) in a frame of refer-noe 

moving at the local (low velocity (compare the definition (25) of uy). The kinetic-energy 
part of the excess energy density E is 


< 2 (pf + P*)( u /i + U ai ) 2 > — nPfUjiUfi, 


which includes an extra term 

I >= 

and E is the sum of expressions (28) and (30), giving 


r , U. 

< P,u ft u 3 , > = < P ' Uji lLp' 

Pf k 


E ~ Et ~' or ’ equivalently, either E r = E— or — = — 

r U> u r u ; ' 


(29) 


(30) 


(31) 


'I he quantity Efu, called action density in Hamiltonian mechanics, is identical In both 
rames of reference, and Equations (24) and (27) tell us that its flux along a ray is constant 

(WF and Lighthill, 1972). 

But Equation (31) shows too that energy is exchanged between (i) the acoustic motions 
relative to the mean flow and (ii) the mean flow itself. For example, where sound waves of 
frequency u> enter a region of opposing flow (or leave a region where the mean flow is along 
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their direction of propagation) the ratio u> r /u> increases and so therefore does E r jE\ the 
sound waves gain energy at the expense of the mean flow. 

The rate of exchange of energy takes the value (1) written down in §1. This is readily 
seen from the laws governing motion in an accelerating frame of reference, which feels 

an inertial force = -(mass) x (acceleration of frame). (32) 

If we use at each point of space a local frame of reference moving with velocity u y, then 
fluid in that frame has velocity u si but is subject to an additional force (32); where, per unit 
volume, mass is pj and the frame’s acceleration is 

9 uji . . t du/i . dui; 

Us} ~dx~ glving force ~ PiUs3 ~foT domg work ~ Pi < u *' u *j > (33) 

per unit time on the local relative motions. This rate of energy exchange (33) proves to be 
consistent with the fact that it is the flux, not of E r but of action E t /uj t , that is conserved 
along ray tubes. 

Energy can be extracted from a mean flow, then, not only by turbulence but also by 
sound waves; and, in both cases, the rate of extraction takes the same form (33) in terms of 
perturbation velocities u a ,. It represents the effect (§1) of the 


mean momentum flux pj < u,iu aj > 


(34) 


or Reynolds stress (Reynolds, 1895) with which either the sound waves or the turbulent 
motions act upon the mean flow. For sound waves, by Equation (28) for E r and by the 
substitution 

ct ki k t , 

u , i = — p s — , mean momentum flux = E r —k ; (35) 

pf k kj 

so that the Reynolds stress is a uniaxial stress in the direction of the wavenumber vector 
having magnitude E r . 

[Note: strictly speaking, the complete 


radiation stress E r 


(bh 
\ * 


+ 



for sound waves 


(36) 


includes not only the momentum flux (35) but also the waves’ mean pressure excess 


\_ 

2 



<P\ 


>- 


izli 

2 Pj 


p] 


■>= 


7 ~ 1 


E r 


(37) 


acting equally in all directions (Bretherton and C.arrett, 1968); however (§1) this isotropic 
component produces no energy exchange with solenoidal mean flows.) 
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5-3. Propagation through sheared stratified winds 


The extremely genera] ray-acoustics treatment outlined above for sound propagation 
through fluids in motion has far-reaching applications (in environmental and, also, in engi- 
neering acoustics) which, however, are illustrated below only by cases of propagation through 
parallel flows, with stratification of velocity as well as of temperature (WF and Lighthill, 
1972). The X \ -direction is taken as that of the mean flow velocity ^(#3) which, together 
with the sound speed c(x 3 ), depends only on the coordinate x 3 . Thus, V replaces u/, in the 
general theory while c replaces cj (and, for atmospheric propagation, x 3 is altitude). [Note: 
the analysis sketched here is readily extended to cases of winds veering with altitude, where 
Uf 2 as well as u /1 is nonzero.] 

Either the basic law (23) or its ray-acoustics form (26) provides, in general, “refraction” 
information in the form of three equations for change of wavenumber; while the single, far 
simpler, Equation (24) is a consequence of, but is by no means equivalent to, those three 
By contrast, in the particular case when uj t and cj are independent of Xi and x 2 , Equations 
(24) and, additionally, (26) in the cases j — 1 and 2 give three simple results, 


uj = constant, k } = constant and k 2 = constant along rays, 

that may be shown fully equivalent to the basic law. 

If now we write the wavenumber (a vector normal to crests) as 

(k 1, k 2 , k 3 ) = (k cos V\ * sin ?/>, k cot 0 ), 


(38) 


(39) 


so that k is its constant horizontal resultant, tf) its constant azimuthal angle to the wind 
direction, and 9 its variable angle to the vertical (Figure 6), and use Equation (25) in the 
form 

u; = c(x 3 )k - f V(x 3 )k] = c ( x 3 ) kc . oseefl + K(x 3 )ic cos 0, (40) 

we obtain the important 


extension sin^ = 


to Snell’s Law 


(41) 


U)K 1 - V (x 3 ) COS V' 

from the classical case (V = 0) when the denominator is a constant. This extended law (41) 
tells us how 6 varies with x 3 along any ray - whose path we can then trace, using Equations 
(26) in the form 


^ = c(x 3 ) cos t/> sin 0 + V (x 3 ), = c( 1 3 ) sin r/’ sin 9 , 


dx 3 

It 


(42) 


= r(jr 3 )cos0, 
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by simply integrating dx\jdx 3 and dx 2 /dx 3 with respect to x 3 . 

It follows that a ray tube covers the same horizontal area at each altitude, so that 
conservation of the flux of wave action /? r /u> r along it implies that the vertical component 

( E T /uj r )(dx 3 /dt ) = E t k~ x sin 9 cos 9 (43) 

of wave action flux is constant along rays; from which, with Equation (28), sound amplitudes 
are readily derived. 

Wind shear is able to reproduce all the main types of ray bending (WF) associated 
with temperature stratification, and ofte>. to an enhanced extent. Roughly, the downward 
curvature of near-horizontal rays in (km) -1 comes to 

3 {V"(x3)cos 4- c'(x 3 )} (44) 

where the velocity gradients are in s -1 and the factor 3(km) -1 s outside the braces is an 
approximate reciprocal of the sound speed (WF and Light hi 11 , 1972). 

Cases when (44) is negative: curvature is upward; its magnitude with zero wind is at most 
0.018 (km) -1 (because temperature lapse rate in stable atmospheres cannot exceed 10°C per 
km, giving c' = -0.006s -1 ) but with strong wind shear can take much bigger values for 
upwind propagation (i/> = 7 r). In either case Figure 7(a) shows how the lowest ray emitted 
by a source “lifts off” from the ground, leaving below it a zone of silence (on ray theory - 
actually, a zone where amplitudes decrease exponentially with distance below that ray). 

Cases when (44) is positive: curvature is downward, as found with zero wind in temperature- 
inversion conditions (e.g., over a calm cold lake) and even more with strong wind shear for 
downwind propagation (ip — 0). Figure 7(b) shows how this leads to signal enhancement 
through multiple-path communication. 

In summary, then, the very familiar augmentation of sound levels downwind, and diminu- 
tion upwind, of a source are effects of the wind’s shear (increase with altitude). 


5.4. Wider aspects of parallel-flow acoustics 

The propagation of sound through parallel flows at wavelengths too great for the appli- 
cability of ray acoustics can be analyzed by a second- order ordinary differential equation. 
Thus, a typical Fourier component of the sound pressure field takes the form 


p 4 (i 3 )e’ (l *"~* ,I, " A!jXj) with p-j- 

dx 3 

[1 k\ + kl 
+ [c 2 (u>-I/fc,) J J 


1 


dp. 


[p(u> - Vk\) 2 dx 3 

P, = o. 


( 43 ) 


Equation (45) can be used to improve on ray acoustics 
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(a) near caustics (envelopes of rays) where it allows a uniformly valid representation of 
amplitude in terms of the famous Airy function, giving “beats” between superimposed 
waves on one side of the caustic, and exponential decay on the other (WF); 

(b) at larger wavelengths by abandoning ray theory altogether in favor of extensive nu- 
merical solutions of Equation (45); and 

(c) to obtain waveguide modes for sound propagation in a two-dimensional duct (between 

parallel planes) (Pridmore-Brown, 1958; Mungur and Cladwell, 1969; and Shankar 
1971). 

On the other hand, in the case of a three-dimensional duct carrying parallel flow V(x 2 , x 3 ) 
m the x, -direction, Equation (45) is converted into a partial differential equation (the first 

term being supplemented by another with d/dx 2 replacing d/di 3 , while k 2 is deleted) which 
is used 


(d) to obtain waveguide modes in such ducts; 

(e) in calculations of propagation of sound through the wide jets - modelled as parallel 
flows - typical (§4.1) of modern aero-engines; and, 

(f) with aeroacoustic source terms included, in certain enterprising attempts at modelling 
jet noise generation and emission (Phillips, 1960 and Mani, 1976). 


6. Acoustic Streaming 

61 - Streaming as a result of acoustic attenuat ion 

Sound waves act on the air with a Reynolds stress (34) even when mean flow is absent (so 
that subscript / becomes subscript zero). The j-component of force acting on unit volume 
of air (WF, pp. 337-) is then 


~ -Q - < PoUsiU SJ >: ( 46 ) 

the force generating acoustic streaming (Lighthill, 1978a). 

However, the force (46) could not produce streaming for unattenuated sound waves; 
indeed, their linearized equations can lv used to show that 


if p M =< -CqPq'p* 


\_ 

2 ' 


oPo«*i«a. > then Fj - 


dp 


M 


d 


dx, < dt ip,Uv) > 


(47) 
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which is necessarily zero (as the mean value of the rate change of a bounded quantity). 
Accordingly, the fluid must remain at rest, responding merely by setting up the distribution 
p M of mean pressure whose gradient can balance the force. [Note: actually, on the ray- 
acoustic approximation (28), p M is itself zero, but the above argument does not need to use 
this approximation.] 

Attenuation of sound waves takes place 

(a) in the bulk of the fluid through the action of viscosity, thermal conductivity and lags 
in attaining thermodynamic equilibrium (Ch. 6); and 

(b) near solid walls by viscous attenuation in Stokes boundary layers. 

All these effects produce forces (46) which act to generate acoustic streaming. It is 
important to note, furthermore, that even the forces due solely to viscous attenuation - 
being opposed just by the fluid’s own viscous resistance - generate mean motions which do 
not disappear as the viscosity p tends to zero (Rayleigh, 1896; Nyborg, 1953; Westervelt, 
1953; and Nyborg, 1965). 

6.2. Jets generated by attenuated acoustic beams 

Attenuation of type (a) produces a streaming motion uj : satisfying 

pu.jidu.fj/dx, = Fj - dp/dxj + pS/ 2 u }] . (48) 

Substantial streaming motions can be calculated from this equation only with the left-hand 
side included (Stuart, 1966); although in pre-1966 literature it was misleadingly regarded as 
“a fourth-order term” and so ignored - thus limiting all the theories to uninteresting cases 
when the streaming Reynolds number would be of order 1 or less. 

We can use streaming generated by acoustic beams to illustrate the above principles. If 
acoustic energy is attenuated at a rate ,6 per unit length, then a source at the origin which 
beams power P along the aq-axis transmits a distribution of 

power Pe"^ X| , and therefore energy per unit length c~* Pc ~^ Zx ; (49) 

which is necessarily the integral of energy density, and so also of the uniaxial Reynold stress 
(35), over the beam’s cross-section. It follows by differentiation that the force per unit 
volume (46), integrated over a cross-section, produces (Lighthill, 1978a) 

a force c~ x P(ie~^ x ' per unit length in the X) -direction. (50) 

At high ultrasonic frequencies the force distribution (50) is rather concentrated, the distance 
of its center of application from the origin being just ()~ x (which at 1 MHz, for example, is 


19 


24mm i„ air). Effectively ( WF, p. .145), the beam applies at this center a total force c- P 
(integral of the distribution (50)). 

The type (WF and Liglithill, 1978a) of streaming motion generated by this concern rated 
force c- ' P depends critically on the value of pc- a sort of Reynolds number squared, 

winch is about 10' P in atmosph. ric air (with P in watts). Streaming of the low-Reynolds- 
number “stokeslet" type predicted (for a concentrated force) by Equation (48) with the left- 
band side suppressed is a good approximation only for P < Kr'HV. For a source of power 
10- H', by contrast, the force c"P generates quite a narrow laminar jet with momentum 
transport r- P, and at powers exceeding 4 X 10- W this jet has become turbulent, spreading 
conically with semi-angle about 15“ and continuing to transport momentum at the rate ' P 

•Such turbulent jets generated by sound are strikingly reciprocal to a classical aeroacoustic 
theme! 

At lower frequencies an acoustic beam of substantial power delivers a turbulent jet with 
a somewhat more variable angle of spread but one which 

at each point X\ carries momentum transport c~ } P(\ — ^]) 

generated by the total force (50) acting up to that point. This momentum transport in the 
jet represents the source’s original rate of momentum delivery minus the acoustic beam’s 
own remaining momentum transport (49). In summary, as acoustic power is dissipated into 
ieat, the associated acoustic momentum transport is converted into a mean motion (which, 
at higher Reynolds numbers, is turbulent: Liglithill, 1978a). 

6 3 - ^reaming around bodi es generated by boundary-layer attenuation 

Sound waves of frequency u, well below high ultrasonic frequencies have their attenuation 
concentrated, if solid bodies are present, in thin Stokes boundary layers attached to each 
body ( WF). Then the streaming generated near a particular point on a body surface is rather 
simply expressed by using local coordinates with that point as origin, with the r-axis normal 
to the body and the .r-axis in the direction cf the invisrid flow just outside the boundary 
layer: the exterior flow. The Stokes boundary layer for an exterior flow 

(f/(j,.y), l /(r,y))e ,w has interior flow 

(C(j \ y), V{x,y))c iut [l - ^ 

[Note that my choice of coordinates makes T(0,0) = 0, and that the expressions (52) become 

identical outside the layer.) The streaming (WF and Ligl, thill, 1978a) is calculated from the 
equation 

i-r - rr + ,*.«/*• - o, (IB) 
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with certain differences from Equation (48) explained as follows: 

(a) the first term is the force (46) generating streaming within the boundary layer; 

(b) we are free, however, to subtract the second, since (see §6.1) it can produce no stream- 
ing, and conveniently, the difference is zero outside the layer; 

(c) gradients in the 2-direction are so steep that the third term dominates the viscous force 
- and, indeed, in such a boundary layer, dominates also the left-hand side of Equation 
(48). 


1 he solution of Equation (53) which vanishes at z = 0 and has zero gradient at the edge 
of the layer is obtained by two integrations, and its exterior value is 


!M) 

where integration extends in practice, not to “infinity,” but to the edge of the layer within 
which the integrand is nonzero. Expression (54) for the exterior streaming is yet again (see 
§6.1) independent of the viscosity p since Equation (52) makes zdz of order p/puj\ and it is 
easily evaluated. 

At x = y = 0 (in the coordinates specified earlier) the exterior streaming (54) has 


4 lr m//dx + 2dV/dy , , r dV/dx 

x — component — U and y — component — U — , (55) 


4u> 


4u; 


with zero 2-component. This is a generalized form of the century-old Rayleigh law of stream- 
ing (which covers cases when V is identically zero). 

For the complete streaming pattern, expressions (55) are, effectively, boundary values 
for its tangential component at the body surface (because the Stokes boundary layer is so 
thin). 1 herefore, any simple solver for the steady-flow Navier-Stokes equations with specified 
tangential velocities on the boundary allows the pattern to be determined. Important note: 
here, the inertial terms in the Navier-Stokes equations must not be neglected, unless the 
Reynolds number R , based on the streaming velocity (55) be of order 1 or less; when, 
however, the corresponding streaming motions would (as in §6.2) be uninterestingly small. 

In the other extreme case when R, is rather large (at least 10 3 ) the streaming motion 
remains quite close to the body (Stuart, 1966) within a steady boundary layer whose di- 
mension (relative to that of the body) is of order R' l/i . This layer is by no means as thin 
as the Stokes boundary layer, but it does confine very considerably the acoustic streaming 
motion. Equations (55) direct this motion towards one of the exterior flow's stagnation 
points, whence the steady-boundary-layer (low emerges as a jet (Figure 8) yet another jet 
generated by sound (Riley, 1987 and Amin and Riley, 1990). 
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Figure 1 

(a) When a source of sound of period T ap- (b) Case of observer at angle 0 to source’s 

proaches an observer at velocity u>, succes- direction of motion at speed V ; then veloc- 

sive wave crests are emitted towards him ity of approach is w = V'cos 0. j 

with separation A = cT — wT. i 



Figure 2. Sounds arriving simultaneously from a source’s far side F and near side N took 
times t and t — r (say) to reach the observer, while the distance FN increased from ( to 
( + wr during the time r between emissions from F and N. 
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Figure 4. A uniformly valid Doppler-effect, approximation. 

Diagram (a) Spare-time diagram for unconverted "eddies” of correlation length ( and 
duration 8. 

Diagram (b) Case of “eddies” converted towards observer at velocity w; being Diagram 
(a) sheared by a distance tv per unit time. Here, lines sloping by a distance r per unit time 
represent emissions received simultaneously by observer. 

Case (i): wfc small. Case (ii): t v/c - 1. Case (iii): intermediate value of ic/r. 

Diagram (c) Average modification factor (Id). 

Acoustic efficiency, obtained by applying this factor to a low-Mach-number 

“quadrupole” efficiency of (say) 10 _:, (V / /r) s 


iSfc 


Figure 5. Equation (17) tells us that any value of uj propagates unchanged along a char- 
acteristic dx\/dt\ = uj. This implies (WF, p. 151) waveform shearing at a uniform rate 
(Diagram (a)). According to nonlinear acoustics (WF, pp. 170-) area-conserving discontinu- 
ities, shown as broken lines in Diagram (a), have to be incorporated wherever necessary tu 
keep the solution one-valued; leading, in the case illustrated, to the N - wave form (Diagram 
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Figure 7. Effects of ray curvature (44) on propagation from a source on horizontal ground. 

Diagram (a) Rays of given upward curvature (due to temperature lapse or upwind prop 
agation) can leave a Zone of Silence (ZS) below the ray emitted horizontally. 

Diagram (b) Rays of given downward curvature' (due to temperature inversion or down- 
wind propagation) can enhance received signals through multiple- path communication. 



Incident 
sound waves 



Figure 8. Plane sound waves incident on a sphere may generate a steady streaming motion 
(Amm and Riley, 1990), concentrated in a relatively thick boundary layer, and directed 
towards one of the exterior flow’s stagnation points - whence it emerges as a jet. 
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